The Astrophysical Journal 

Preprint typeset using L^T^X style cmulatcapj v. 10/09/06 



ON THE INTERPRETATION OF THE AGE DISTRIBUTION OF STAR CLUSTERS IN THE SMALL 

MAGELLANIC CLOUD 

Mark Gieles 

Astronomical Institute, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands and 
European Southern Observatory, Casilla 19001, Santiago 19, Chile 



Henny J.G.L.M. Lamers 

Astronomical Institute, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands 

Simon F. Portegies Zwart 

Astronomical Institute 'Anton Pannekoek', University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands, Section 
Computational Science, University of Amsterdam, Kruislaan 403, 1098 SJ, Amsterdam, The Netherlands 

Received 2007 May 9; accepted 2007 June 8 



ABSTRACT 

We re-analyze the age distribution (dN/dt) of star clusters in the Small Magellanic Cloud (SMC) 
using age determinations based on the Magellanic Cloud Photometric Survey. For ages younger than 
3 x 10 9 yr the dN/dt distribution can be approximated by a power-law distribution, dN/dt oc t^ 13 , with 
—j3 — — 0.70±0.05 or —j3 = — 0.84±0.04, depending on the model used to derive the ages. Predictions 
for a cluster population without dissolution limited by a U-band detection result in a power-law dN/ dt 
distribution with an index of ~ —0.7. This is because the limiting cluster mass increases with age, 
due to evolutionary fading of clusters, reducing the number of observed clusters at old ages. When a 
mass cut well above the limiting cluster mass is applied, the dN/dt distribution is flat up to 1 Gyr. 
We conclude that cluster dissolution is of small importance in shaping the dN/dt distribution and 
incompleteness causes dN/dt to decline. The reason that no (mass independent) infant mortality of 
star clusters around ~ 10 — 20 Myr is found is explained by a detection bias towards clusters without 
nebular emission, i.e. cluster that have survived the infant mortality phase. The reason we find no 
evidence for tidal (mass dependent) cluster dissolution in the first Gyr is explained by the weak tidal 
field of the SMC. Our results are in sharp contrast to the interpretation of Chandar et al. (2006), who 
interpret the declining dN/dt distribution as rapid cluster dissolution. This is due to their erroneous 
assumption that the sample is limited by cluster mass, rather than luminosity. 

Subject headings: galaxies: individual (Small Magellanic Cloud) — galaxies: star clusters — stars: 
formation 



1. INTRODUCTION 

Star clusters are often assumed to be tracers of the star 
formation history of their host galaxy. Recent studies on 
the star formation rate in the solar neighborhood have 
revealed th at the majority of stars form in a cluste red en- 
vironment (|Clarke et alJl2000t lLada fc Ladall2003n . while 
only a few percent of stars in the solar neighborhood are 
in clusters. Understanding the process of cluster dissolu- 
tion is, therefore, of key importance if one wants to make 
a meaningful translation from the observed age distribu- 
tion of star clusters to a star formation history of their 
host galaxy. 

Recent theoretical work suggest that a large fraction 
(50-90%) of the star clusters disperse a few Myrs af- 
ter formation due to the expulsion of residual gas by 
the stellar winds and supernovae of m assive stars (e.g. 
iGoodwinl [19971 : iBoilv fc Kroup"al [2001 . If the star for- 
mation efficiency is independent of cluster mass (M c \), 
then this "infant mortality " is also independent of M c \ 
(jGoodwin fc Bastiad[2006li . 

The observed age distribution of star clusters in M51 
(jBastian et al.l 120051 ) indicates that ~ 70% of the clus- 
ters disperse on a time-scale of ~ 10 Myr after formation, 
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roughly independent of M c \. iFall et a l. (2005) claim that 
in the "Antennae" galaxies this process removes roughly 
90% of the clusters each age dex during the first Gyr. 
This time-scale is far too long to be explained by the gas 
expulsion scenario. 

The M c \ independent dissolution due to gas expulsion 
is in sharp contrast to diss olution due to two-body re- 
laxation in a tidal field (e . g. IChernoff fc Wcinbcr gjl!99 
Vesperini fc Heggiel Il997t [T akahash i fc Portegies Zwar 



2000: iBaumgardt fc M akino I2003D, or external 
perturbations by the disk ( Ostriker et al.l 119721 : 



iGnedin fc Ostrikeil Il997f) or gian t molecular clouds 
(jSpitzerl 119581 : iGieles et al.l l2006b[ ). which all have a 
more destructive effect on low mass clusters than on 
high mass clusters. 

The observed age distributions of star clusters can be 
used to disentangle the two aforementioned disruption 
processes: mass independent or mass dependent cluster 
dissolution. If the disruption time (t<jis) has a power-law 
dependence on the cluster mass, then the power-law in- 
dex can be derived from the slope of the age distribution. 
This holds for mass limited and magnitude limited clus- 
ter samples, since at old ages disru ption will usually dom- 
inate over evolutionary fading (|Boutloukos fc Lamersl 
120031) (BL03). On the other hand, in the case of a M cl 
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independent infant mortality process the age distribution 
of a mass limited cluster sample results in an age distribu- 
tion of the form dN/dt oc i -1 , where the index of —1 ap- 
plies to t he case where 90% of the clusters dissol ves each 
age dex (jFall et all 120051 : IWhitmore et~aTll2007f ). If the 
sample is magnitude limited, then the age distribution 
will be steeper. This is because fading and disruption 
both remove a certain fraction from the dN/ dt distribu- 
tion. 

Recently, IChandar et ail (j2006ft (CFW06) studied 
the age distribution of clusters in the Small Magel- 
lanic Cloud (SMC) based on cluster ages derived by 
iRafelski fe Zaritskvi (|2005ft (RZ05). They claim that the 
age distribution is consistent with a power-law with in- 
dex — 1 for t < 3 Gyr and that the sample is "reasonably 
complete down to 10 3 M Q over this age range" . They 
conclude that their observations are consistent with in- 
fant mortality being at work for ~ 1 Gyr and independent 

Of M c i. _ 

In this paper we address the interpretation of the 
dN/dt distribution of star clusters in the SMC and the 
pitfall caused by detection incompleteness that can lead 
to the fallacious conclusion that mass independent clus- 
ter disruption shapes the dN/dt distribution. 

The structure of the paper is as follows. In § 2 we dis- 
cuss the observations from which the SMC cluster sam- 
ple was derived and the mass and age determinations of 
the clusters. Because the interpretation of the age dis- 
tribution of the cluster sample depends heavily on the 
(in)completeness of the sample, we investigate different 
incompleteness scenarios in § [3] In § 0] we compare the 
resulting age distribution to predictions and argue that 
the shape of the age distribution it can be explained by 
evolutionary fading, without involving mass independent 
dissolution. In §[S]we compare our results to other stud- 
ies and our conclusions are outlined in § [6] 

2. OBSERVATIONS 

The clusters analy zed by RZ05 were identified by 
iHill fc Zaritskvi (|2006l ) (HZ06). Photometric data in the 
UBVI bands fr om the Magellanic C louds Photometric 
Survey (MCPS; IZaritskv et ai1ll997f ) were used to find 
clusters. Stellar density images based on the photomet- 
ric catalog of SMC stars were constructed by counting 
the number of stars with V < 20.5 in squares of 10". We 
note that for the photometry all stars in the catalogue 
were used, but only for clusters that were identified in 
the sample limited to stars with V < 20.5. 

RZ05 derived ages by comparing the U — B, B — V 
and V - I colors from clusters using S TARBURST99 
(hereafter SB 9 9, iLeitherer et all 1 1999ft and GALEV 
(jSchulz et all I2002T lAnders fc Fritze-v. Alvenslebenl 
12003ft models. We adopt their ages, which have been 
derived from the photometry with the mean foreground 
extinction correction for the SMC of E(B — V) =0.09, 
but without a correction for local extinction. RZ05 
tried to include extinction as a free parameter in the 
age fitting method, but conclude that the scatter in the 
photometry is too large to improve the age estimates by 
including extinction as an additional parameter. They 
take extinction into account in an analytical way 1 , but 
again conclude that the spread of observed clusters 

1 RZ05 adopted the mean extinction measured by 
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Fig. 1. — Age-mass diagram of 195 clusters in the SMC based 
on ages derived by RZ05 with the GALEV models (bullets) and 
the SB99 models (squares). The 80% limit as a function of age 
[M^'° (t)) is derived from the data (see text for details) and shown 

as triangles for the SB99. A polynomial fit to M^^°(t) is shown 
as a full line. The vertical dashed lines indicate the ages where the 
M* 0% (t) curve intersects with logM c i = 3 and logMd = 3.5. 

colours around the models does not reduce. Since this 
artificially imposed extinction as a function of age could 
introduce systematic trends in the dN/dt distribution, 
combined with the conclusion of RZ05 that this method 
did not reduce the scatter of the cluster colors, we prefer 
the first order age estimates of RZ05. 

Adopting the ages of the clusters based on their pho- 
tometry, we derived initial masses of the clusters, i.e. 
corrected for mass loss due to stellar evolution, indepen- 
dently using SB99 and GALEV models with Z = 0.004. 
For the SB99 models we adopt identical settings as 
RZ05, i.e. standard mass loss, the full isochrone mass in- 
terpolation, a Salpeter initial mass function from 0.1 M 
to 100 M Q . 

In Fig. [T] we show the resulting age-mass diagram for 
the 195 clusters for which an age estimate is available. 
For both samples an increase in the upper and lower 
mass boundaries with age can be seen. The increasing 
upper mass boundary is due to a size of sample effect, 
since the x-axis is the logarithm of age. For a constant 
cluster formation rate, more clusters are formed in each 
dex for increasing logi. For a power-law cluster initial 
mass function (CIMF) with index —2 the upper bound- 
ary is expected to increase as a straight lin e with slope 
+1 (jHunter et alj|2003t iGTeles et alj|2006aj) . m reason- 
able agreement with Fig. [1] 

The sloping lower mass boundary is due to incomplete- 
ness, most likely because of evolutionary fading of clus- 
ters with age. Figure [1] clearly shows that the limiting 
cluster mass in the sample of RZ05 is not constant, con- 
trary to the assumption made by CFW06 in their anal- 
ysis of the age distribution. Note that CFW06 do not 
present an age-mass diagram in their work to support 

I Harris fc Zaritsky (2004|) for stars younger than 10 Myr and 
older than 1 Gyr. For intermediate ages they interpolated the 
extinction value. RZ05 then corrected the colors for the clusters 
for which an age estimate is already available, and recalculated 
the best fit with the predicted clusted models. 
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their assumption, but rely on the assumption that the 
sample "is likely to be approximately surface brightness 
limited, rather than luminosity limited" . In § [3] we in- 
vestigate the possible causes for incompleteness and the 
effect on the dN/ dt distribution in more detail. 

Because of the method used by HZ06 to define the 
cluster sample, the sample of RZ05 is not expected to 
be magnitude limited only (see § [3|. In an attempt to 
quantify the increasing lower M c \ as a function of age, we 
empirically derive the mass limit above which 80% of the 
clusters are more massive (M 80% (t)). The limit is chosen 
at 80% since we find this to be the point where the mass 
function (at different ages) turns-over. This suggests the 
sample is relatively complete above this limit and highly 
incomplete below it. 

At intervals of 0.25 dex in logi we count the number 
of clusters in a bin with a width of 0.5 dex. The masses 
in the bin are sorted and we count from the highest mass 
until we have reached 80% of the total number in the 
bin. Note that this M^ % {t) is not the same as an 80% 
completeness limit, or the detection limit, as usually de- 
rived from artificial cluster experiments, since we here 
have no information about clusters that did not make it 
into the sample. However, if we assume that the shape of 
the completeness curve (i.e. the function that describes 
the fraction of clusters that is retrieved as a function of 
mass), is not dependent on age, we can safely assume 
that the evolution of with age is the same shape 

as that for the completeness limit. 

The result of M®® % (t) is shown as triangles in Fig. Q] A 
polynomial fit is shown as a full line. Note that the curve 
is located relatively high above the lower limit. If the 
sample had a "hard" detection limit and the cluster IMF 
had an index of —2 and the cluster dissolution would be 
mass independent, we would expect that at each age the 
MS (t) line is only a factor 1.25 higher in mass than the 
lower limit. Fig.[T]shows that the 80% limit is at about an 
order of magnitude higher mass than the absolute lower 
limit. This implies the data is highly incomplete for low 
mass clusters. In the next section we will compare the 
M^ % {t) curve to different predictions. 

3. INCOMPLETENESS: LIMITED BY MASS, BY 
LUMINOSITY OR A COMBINATION OF BOTH? 

The interpretation of the empirical age distribution of 
the SMC clusters depends strongly on how incomplete- 
ness affects the sample. Therefore, we compare the em- 
pirically derived (t) curve to what is expected from 
the selection procedure of HZ06 and RZ05. 

3.1. A sample limited by the number of stars in a 
cluster 

From the stellar catalogue HZ06 selected stars with 
V < 20.5 mag. They constructed number density images 
from which over-densities were detected. This implies 
that only clusters with enough stars with V < 20.5 could 
end up in the catalogue. 

We first assume that the number of stars with V < 20.5 
(N(V < 20.5)) determines the limiting cluster mass. We 
use the evo lutionary isochrones of the Padova models for 
Z = 0.004 (jBertelli et al.l[i~99i IChardi et al]ll99dl20"0Ch 
which wer e converted to the U BVRIJHK photometric 
system by iGirardi et all (|2002t ). We assume a Salpeter 
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Fig. 2.— The empirically derived M^ 0% (t) curve above which 
80% of the clusters are located as derived in §[2]and Fig.[T](full line). 
Also shown are the two predicted lines for the cluster models with 
a limiting mass evolution based on the assumption of a constant 
number of stars brighter than V < 20.5 (dotted line) and on the 
assumption of a limiting cluster magnitude at My < —4.5 (dashed 
line). 

stellar initial mass function (IMF) with lower mass of 0.5 
Mq. Stars below this mass do not reach V < 20.5 within 
the Hubble time. The number of stars with V < 20.5 is 
counted for cluster models of different ages and masses 
and we select only clusters with N(V < 20.5) greater 
than some fixed number. We adopted N(V < 20.5) = 25 
because the resulting mass-age relation agrees more or 
less with the location of the M c 8 j 0% (t) line in Fig. 2. 

In Fig. [5] we show the result as a dotted line. The curve 
is slowly increasing for logt < 9 and then rises. This 
transition occurs at the age where the main sequence 
turn-off drops below V — 20.5. Below that age a clus- 
ter of a given mass has nearly constant N(V < 20.5), 
since the brightness of stars on the main sequence is not 
varying much and the fraction of stars in N(V < 20.5) 
that is on the main sequence is high. When the turn- 
off drops below V = 20.5, N(V < 20.5) consists of RGB 
and AGB stars mainly. Because the fraction of RGB and 
AGB stars in a cluster is small, but is responsible for the 
majority of the cluster luminosity, the minimum mass of 
a cluster with N(V < 20.5) = 25 steeply increases with 
age. 

3.2. A sample limited by the total cluster luminosity 

After HZ06 selected their clusters from the stellar den- 
sity image, as described in § [2j a visual inspection on the 
original images was done. Only the unambiguous clus- 
ters were kept. After this, they fitted surface brightness 
profiles to the remaining sources. Only when a reliable 
surface brightness fit and photometry could be extracted, 
the source was kept in their sample and called a "star 
cluster" . The total number of sources found from the 
stellar density images was a factor of 4 higher than the 
final number of star clusters in the catalogue (Zaritsky, 
private communication). This indicates that accurate 
photometry was a strongly limiting factor in the selec- 
tion of the clusters. Since this last selection step relies 
on the total luminosity of a cluster we may assume that 
the sample is mainly limited by the magnitude of the 
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clusters. To test this assumption, we predict what is ex- 
pected for the limiting mass as a function of age for a 
constant magnitude limit. 

Star clusters are formed on a time-scale short enough 
that their photometric evolution can be well described 
by Simple Stellar Population (SSP) models. From these 
models it follows that the flux of a star cluster scales ap- 
proximately as a power-law with age (F\ oc t~^), with 
£ positive. The limiting cluster mass as a function of 
age (M^™(t)) of a sample that is limited by an absolute 
magnitude at a certain wavelength (M} lm ) can be ex- 
pressed in the magnitude evolution from an SSP model 
(Mf p (i)): 

logAf'j m (t) = logAfff - 0.4 [i\4 lm - Aff sp (t)] , (1) 

where is the initial mass of the cluster that is de- 

scribed by the SSP model. For the GALEV and SB99 
version that we use M^ff = 10 6 M Q . When taking M^ s f 

as a constant M^ m (t) represents an initial mass, i.e. be- 
fore mass loss due to stellar evolution. 

Since logM^ p (t) and M^ m from Eq. [Hare constant 

and 0.4M| SP oc - log Fx = Clogt, we can rewrite Eq. Q] 

as 

log M^(t) = (logt + C, (2) 

with C = log Af c S jSP _ 0.4M| im (BL03). 

In Fig.[2]we show the result of Eq.[T]as a dashed line for 
the SB99 models with an adopted limiting magnitude of 
My m = —4.5. This value was chosen because the result- 
ing limiting age-mass relation agrees more or less with 
the location of the empirically derived Mjj (i) limit of 
Figs. Q] and [2] Notice that this prediction matches the 
empirically derived Af' I 1 m (t) quite well. We especially 
note that this predicted relation fits the empirical one 
much better than the N(V < 20.5) curve. 

We conclude that the cluster sample of RZ05, used to 
derive the age distribution of SMC clusters, is mainly 
magnitude limited. 

4. THE AGE DISTRIBUTION OF SMC CLUSTERS 
4.1. The observed age distribution 

The uncertainties in the age determinations of RZ05 
are quite large, mainly due to the photometric uncer- 
tainty of faint clusters. We take these uncertainties into 
account when constructing the dN/dt distribution. 

We represent the contribution by each cluster to the 
age distribution as an asymmetric Gaussian, with the 
lower width (<7_) and the upper width (<r+) correspond- 
ing to the minimum and maximum age (t_ and t+) de- 
rived by RZ05. RZ05 based these latter values on the 
90% confidence range from their x 2 results of the pho- 
tometric age determination. This corresponds to 1.6er. 
We assume a Gaussian spread in \ogt, since the errors 
in the age determination are roughly constant in logt 
(mean A logt ~ 0.25). The definition of a± is then 
a± = | logt± — log 1 1/1. 6. This results in an asymmetric 
Gaussian profile for each cluster, where the left and right 
side both have a surface of 0.5, such that the total contri- 
bution of each cluster to the age distribution is 1 . Some 
ages derived with SB99 have cr_ — or a + = 0. In that 



case we adopted a minimum uncertainty of a± = 0.05. 
When all clusters are added in a large log t array of equal 
intervals in logt, each value of dN (i.e. the number of 
clusters in that age interval) is divided by the width of 
the bin in linear age, to construct the dN/dt distribution. 

In the left and right panel of Fig. [3] we present the 
smoothed dN/dt result for the GALEV and SB99 mod- 
elling, respectively. The full line represents the full data 
set. We estimated Poisson errors by counting the number 
of clusters in bins of width 0.25 dex, corresponding to the 
mean uncertainty in the log age values of RZ05. The la 
Poisson errors are shown as a grey shaded region. Note 
that we do not normalize the cluster age distribution to 
the stellar age distribution, as was done by RZ05. This 
to compare our results to the results of CFW06, who also 
construct the unnormalized cluster age distribution. 

This representation of the dN/dt distribution takes 
into account the increasing age interval for increasing 
width of the age intervals. The number of clusters per 
logarithmic age bin actually slightly increases with logt, 
which reflects the small increase in number density that 
can be seen in Fig. [T] Note that dN/dt oc t~ l implies a 
constant number of clusters in constant logt bins. 

In Fig. [3] we also show a binned dN/dt histogram as 
filled circled with small error bars. The sizes of the bins 
are chosen such that there are 17 clusters in each bin (fol- 
lowing |Mafz5ciMnii5!Tjbeda 2005). This results in a 
total of 10 bins over a range of ~ 4.5 age dex. The ob- 
tained average bin width is then 0.35 dex, which is larger 
than the mean error in the age determination making it 
suitable for a fit where we only take into account the un- 
certainties in the direction of the y-axis. The first bin 
starts at the age of the youngest cluster and the last bin 
ends at the age of the oldest cluster. This way of binning 
the data results in a dN/dt distribution very similar to 
the Gaussian smoothed version (Fig. [3]). 

We fit a straight line to log(dAT/df) vs. logt for 
t < 3 Gyr using a x 2 -error statistics minimization with 
bin weights (Wj) depending on the standard deviation of 
each bin (<jj) as Wi = 1/of , which for this representation 
is Wi — In 2 (10) Ni. We find a slope, i.e. the power-law 
index, of -0.70 ± 0.05 if we use the GALEV results and 
-0.84 ± 0.04 if we use the SB99 result. The mean index 
is thus 0.77 ±0.07. 2 

We conclude that the observed age distribution in the 
age range of 7 < logt < 9.5 can be approximated by a 
power-law with index —0.77 ± 0.07. 

4.2. The predicted age distribution 

Assume that a cluster population formed at a constant 
rate and with a power-law CIMF with index —a. Then, if 
there is no dissolution, the age distribution of all clusters 
can be acquired by integrating over all masses from M^.\ m 
(from Eqs. Q] or |2J| to Af™ ax (BL03): 



r lvl c\ 

dN/dt = / SM c r a dM ci 

2 If wc ignore the bin weights in the fit, we find uncertainties on 
the indices of 0.20 and 0.08 for the fits on GALEV and SB99, re- 
spectively. This could be the explanation for the large uncertainty 
of 0.15 found by CFW06 in their index of 0.85. 
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Fig. 3. — Top: The age or dN/dt distribution from GALEV (left) and SB99 (right) results. Full line: each cluster was represented by a 
Gaussian profile with cr± corresponding to the minimum and maximum uncertainty in the age determination. The grey area corresponds 
to the la uncertainty. Dots with error bars: the age distribution shown as a histogram where the horizontal lines indicating the variable 
width of the age bins, where each bin contains 17 clusters (see text for details). Thick dashed line: the predicted dN/dt distributions for 
a sample limited by M^^° (t) of Fig. [T] The thick dotted line represents a cluster sample limited by M s ®^° (t) combined with the mass 
independent disruption model which was suggested by CFW06. This distribution declines much more rapidly than the data. Bottom panel: 
the ratio between the observed dN/ dt distribution and the predicted one for the sample limited by M^^° (t) . 



s 



1 - a 

rlim 



[( M max)l-a _ (M'i" 1 ) 1 - 



(3) 

In the last 
M'j m . If the 



logM cl 



ocl/M£ 

where S describes the cluster formation rate 
steps we have used a — 2 and M™ ax >> 
sample is magnitude limited, as we have shown in § O 
with M^\ m oc v> then the age distribution scales with t as 
dN/dt oc t~^. So we expect a slope — ( in the logarith- 
mic representation of the dN/dt distribution (Fig. [3]). On 
the other hand, for a sample that is mass limited (cor- 
responding to Af'j 111 being constant) without dissolution 
the age distribution would be flat. 
The value of C can be determined by approximating 
„ S j SP (i) vs. \ogt following from SSP models by a 
straight line. For ages smaller than 3 Gyr, the value of 
C derived from the GALEV models are [1.02, 0.89, 0.72, 
0.63, 0.58] for U, B, V, R and /, respectively For the 
SB99 models we find [0.93, 0.78, 0.66, 0.60, 0.54] for the 
same filters, respectively. So, for a cluster sample formed 
with constant formation rate, not affected by disruption 
and limited by a detection in the V band, the predicted 
age distribution is a power-law with index ~ —0.7, in 
agreement with the observations (§ 14. ip . 

Using the empirically derived Ml\ m (t) from Figs. [1] and 
[51 we predict the dN/dt distribution for a cluster sample 
that is not affected by any disruption and thus declines 
as [M^it)}- 1 (thick dashed lines in Fig. [3]). The predic- 
tion and the observed data are normalized at logt = 8.5. 
The prediction describes the overall decline and the de- 
tails in the shape of the observed dN/ dt distribution very 
well. This suggests that there are no signs for cluster dis- 
ruption present in this data set and the declining dN/dt 
distribution can be explained by incompleteness only. 

To illustrate the combined effect of disruption and in- 
completeness, we also show the dN/dt distribution for 
a sample where we apply the 90% reduction of clusters 
each age dex, as suggested by CFW06 (thick dotted line) , 



also normalized to the data at logt = 8.5. We stress that 
the mass independent disruption model of CFW06 works 
in addition to incompleteness, such that both make the 
dN/dt distribution steeper. Therefore, a fallacious as- 
sumption that a cluster sample is mass limited, while it 
is in reality luminosity limited, will lead to a confusion 
between disruption and incompleteness. 

The nice agreement between the observed dN/dt dis- 
tribution and the [M®® % (i)] _1 curve implies that the ob- 
served age distribution of the SMC clusters with ages 
logt/yr < 9.5 can be explained by evolutionary fading 
only, without the need for mass independent dissolution. 
The fact that there is no evidence for infant mortality 
at young ages (< 30Myr) is most likely due to a bias 
of clusters without nebular emission (HZ06), which typi- 
cally have such young ages. We expect that the effect of 
gas removal must destroy young (few Myrs) clusters in 
the SMC, but we conclude that the sample under discus- 
sion here is not suitable to investigate this phenomenon 
in detail. In a future study we will investigate the infant 
mortality rate of yo ung star clusters i n the SMC using 
different techniques (jGieles et al.ll2007t) . 

Also tidal evaporation does not affect the dN/dt dis- 
tribution for SMC clusters with ages < 3 Gyr. This 
agrees with the long survival times (> 6 Gyr for masses > 
10 4 M©) predicted by iV-body sim ulations of clusters dis- 
solvin g in weak tidal fields (e.g. iBaumgardt fc Makinol 
l200l . 

4.3. The age distribution of a mass limited sub-sample 

The results of CFW06 are based on the assumption 
that the RZ05 sample of clusters is mass limited. To 
show the effect of a mass limit on the age distribution we 
compare the cumulative age distribution of the clusters 
with and without mass limit to the predictions. Such a 
cumulative distribution is preferred over constructing a 
histogram when only a small number of clusters is avail- 
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able, as is the case when applying a (relatively high) mass 
cut (see Fig. [T]). 

Fig. 2] shows the cumulative age distributions divided 
by the age range from the minimum age in the sample 
(i m in): N(t' < t)/(t — frmin). This representation has the 
advantage that it resembles a dN/dt distribution. The 
full SMC cluster sample is shown as thick full lines and 
the sub-sample of clusters with masses above 3 x 10 3 M© 
is shown as thick dashed lines, using the masses and ages 
derived by using the GALEV models (left panel) and 
the SB99 models (right panel). The grey areas of each 
distribution indicates the la uncertainty due to Pois- 
son statistics. Note that the sub-sample with the mass 
cut is flat up to ~ 1 Gyr, while the full sample is de- 
clining. Signs of flattening are also present in the dN/dt 
histogram of CFW06 (their Fig. 1) for a sub-sample with 
M > 1O 3 M0. The authors also note this flattening and 
argue that "the mass-limited sample is somewhat shal- 
lower than the one constructed from the entire sample, 
since there are relatively few clusters more massive than 
10 3 Mq at very young ages." While this is indeed the 
correct explanation for why the dN/dt distribution of a 
mass-limited sub-sample is flatter, it should not be the 
case in their suggested scenario. In their scenario, the 
number of clusters in log t bins of equal width should be 
constant. Using the GALEV results, we find 36 clusters 
with masses above 3 x 10 3 M Q , and ages between 10 8 75 yr 
and 10 9 25 yr, while there are clusters above this mass 
limit with ages between 10 6 75 yr and 10 7 25 yr. Accord- 
ing to the model of CFW06 these numbers should be the 
same. In addition, we show in Fig. [1] that a mass cut at 
10 3 M Q is not high enough to be safely above the increas- 
ing limiting mass up to Gyr. A mass cut at 3 10 3 M Q 
would be a safer choice. Therefore, we will consider this 
limit. 

We compare these cumulative distributions with those 
predicted for the power-law fits to dN/dt. If dN/dt — 
Ct~ v , then the N(t' < t) depends on t as 



N(f <t) = 



C 



1 



V 



1-17 



- t 



1-17 



if V ± 1. (4) 



We adopt i min = 3Myr, since this is the youngest age 
present in the SSP models. 

In § l4.1l we have shown that for the full sample r\ — 0.70 
and j] = 0.84 for GALEV and SB99, respectively. The 
resulting predicted cumulative distribution (Eq.[?]), again 
normalized to (i— i m in), is shown by the thin full line. No- 
tice that it agrees with the empirical distribution within 
about 2cr. For a sample that is mass limited above the 
detection limit and does not loose clusters by dissolution, 
i.e. rj = 0, the expected distribution will be flat in the 
representation of Fig. 0J The thin dotted lines in Fig. 
0] show the predicted distributions. Again, these agree 
with the observed cumulative distributions up to an age 
of ~ 1 Gyr, corresponding to the maximum age for which 
the sample of cluster with log M = 3.5 is not affected by 
the M*° % (t) line (Fig. Q}. Beyond this age the distri- 
bution declines in a similar way as the full distribution, 
which is because even for the relatively high mass cut, 
the sample is affected by incompleteness in the same way 
as the full sample. 

5. COMPARISON TO OTHER STUDIES 



We first compare our results to other studies on the 
SMC star clusters. This is particularly interesting since 
there appears to be a large variation between the different 
data sets and their interpretations. 

One of the first detaile d analy ses of the age distribu- 
tion was done bv IHodgel (|1987| ). The age distribution 
presented in that work looks quite different from the one 
shown here and in RZ05 and CFW06. It is flat up to 
almost 1 Gyr and then has a steep drop. This is quite 
simila r to what is found in a recent study bv lChiosi et al.l 
(2006) (their Fig. 10) who also present an age distribu- 
tion for SMC clusters that is nearly flat the first ~ 10 8 yr. 
Hodgd (119871) dete cts stars down to B ~ 22 - 23 and 
Chiosi et alT(|2006l ) even down to V ~ 24. Both stud- 



ies derive their ages with isochrone fitting techniques. 
This makes it possible to determine the age even if only 
a few stars are available. Assuming that at young ages 
(<few lOOMyrs) all clusters are detected, a flat distri- 
bution is expected until the age that the fading or dis- 
ruption starts to remove clu sters from the sample. Using 
the sample of IHodgel (|1987f) BL03 argued that the disso- 
lution time of clusters in the SMC is very long and of the 
order of 8 Gyr for a cluster of 10 4 M Q . This agrees with 
the fact that we do not find evidence for dissolution up 
to a ges of about 3 Gyr from the RZ05 sample. We note 
that lChiosi et a l. (2006) present 2 age distributions. The 
one we refer to here is for clu sters which were clas sified 
C (^genuin e star clusters) by iBica &; Dutral (|2000D . In 
their Fig. 7 lChiosi et al.l (|2006j ) also present an approxi- 
mately flat age distribution, defined as diV/dlogi, which 
corresponds to dN/dt oc t -1 . However, we note that 
this sample includes all large OB associations, with sizes 
much larger (> 20 pc) than typical star clusters (~ 5pc). 

HZ06 and RZ05 consider only stars with V < 20.5 for 
the selection of their clusters. This cut in the stellar sam- 
ple is ~ 2 — 3 magnitudes brighter than the one on the 
stellar sample mentioned before. In ad dition , the b road 
band photometry of RZ 05 was based on lKingj (|1962f ) and 
EFF (jElson et al.lll987h profile fits to the surface bright- 
ness profile of clusters. More stars are needed for a good 
profile fit to the surface brightness profile than for an 
isochrone fit, which results in les s clusters in the s ample 
of RZ05 than in th e samples of lElson et al.1 (|1987h and 
IChiosi et al.l |2006). This idea is supported by th e larger 
numb er of clusters found in th e deeper studies of [Hodge 
(fl98l and IChiosi et~aTl (p006h . viz. 327 and 311, respec- 
tively. This supports our conclusion that incompleteness 
effects are more important in the sample of RZ05. 

Comparing our results to those of CFW06, based on 
the same cluster sample of RZ05, we note that our in- 
terpretation of the age distribution is drastically differ- 
ent. CFW06 assumed that the cluster sample of RZ05 is 
mass limited. They fit the age distribution of the clus- 
ters with a power-law of index —0.85 ± 0.15, and argue 
that this is consistent within la with mass independent 
dissolution in which the number of clusters decreases by 
a factor 10 per age-dex, up to an age of about 3 Gyr. 
We have shown that the index of the age distribution is 
-0.77 ±0.07 (§|HIJ), i.e. consistent with -1 only within 
3.3<7. Moreover (and more importantly), we have shown 
in § 2 that the cluster sample is not mass limited, as is 
already clear from the distribution of the clusters in the 
age-mass diagram of Fig. [TJ We have given arguments 
that the cluster sample is approximately magnitude lim- 
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10 6 10 7 10 s 10 9 10 10 10 s 10 7 10 s 10 9 10 10 

t [yr] t (yr) 

Fig. 4. — The cumulative age distribution (N(t' < t)) normalized to (t - t min ) for the GALEV (left) and SB99 (right) results. The 
total sample is shown with thick full lines, with the results of the power-law fits overplotted as thin lines. Eq. [4] was used to relate the 
power-law function to N(t' < t). The resulting distributions after a mass cut at log M = 3.5 is applied are shown as thick dashed lines. 
The predictions for a flat dN/dt distribution, up to the age for which the sample is expected not to be affected by incompleteness (~ 1 Gyr, 
see Fig.fTJ, are also flat in this representation (see text and Eq.[4]for details) and are shown as dotted lines. 



ited. This follows both from the way in which RZ05 have 
defined their cluster sample and from the observed 80% 
limit of the stars in Fig. [1] This limit has a very sim- 
ilar slope as that expected for a sample that is limited 
by the visual magnitude (Fig. [2]). The observed index 
of —0.77 ± 0.07 is very similar to the one predicted for 
a sample is magnitude limited in the U-band (0.72 for 
GALEV models and 0.66 for SB99 models) if the num- 
ber of clusters decreases only by evolutionary fading be- 
low the magnitude limit. We thereby rule out the need 
for the infant mortality scenario that proceeds up to 3 
Gyrs as interpreted by CFW06. 

From theoretical arguments it is expected that in- 
fant mortality works on a much shorter time-scale (~ 
10 — 20 Myr) , due to the removal of gas t hat has not been 
used to form stars (e.g. lGoodwmlll997t iBoilv fc Kroupal 
2003). The lack of a strong bump in the age distribu- 
tion of the SMC clusters at young ages (< 30 Myr), as 
was found for the clusters i n M51 dBastian "et~alll2005l) 
and the Antennae galaxies ([Fall et al J 12005( 1 . is proba- 
bly because the SMC sample of HZ06 has a bias towards 
clusters without strong nebular emission. 

6. CONCLUSIONS 

We conclude that the observed age distribution of the 
RZ05 sample of the SMC clusters agrees with that pre- 
dicted for a cluster sample that is magnitude limited and 
in which the age distribution in terms of dN/dt is de- 
creasing with age by evolutionary fading of the clusters. 
There is no need to invoke a mass independent infant 
mortality that extends to 3 Gyr and destroys 90% of the 
clusters every age dex, as suggested by CFW06. In fact, 
for a magnitude limited cluster sample with an extended 
independent mass infant mortality, the age distribution 
would decrease much steeper with age than observed. 

Interestingly, our result also explains the nearly flat age 
distribution of the field stars by RZ05, because not many 
clusters are dissolved between 10 7 and 10 9 yrs, i.e. after 
the infant mortality phase. This is in contrast to the in- 
terpretation by CFW06 of a very long mass independent 
disruption phase, because that model predicts a rapid in- 
crease in the age distribution of the field stars due to the 



continuing contribution of stars from dissolving clusters 
to the field. Such an increasing age distribution of the 
field stars is not observed. 

In this paper we demonstrate that the interpretation 
of the cluster age distribution is very sensitive to in- 
completeness effects as a function of cluster age and to 
how these are taken into account. For extra-galactic 
(slightly resolved) star clusters various tests methods 
have been developed to quantify (in)completeness and 
its dependence on c l uster luminosity a nd radius (see e.g. 
iBastian et aT1l2005t iMora et al J 12007( 1 . For nearby clus- 
ter samples such as the one of the SMC discussed here, it 
should be possible to construct completeness curves as a 
function of structural parameters of the clusters as well. 
By simulating star clusters of different ages and masses, 
thereby taking into account the luminosity evolution of 
the individual stars, and adding them to the catalogues 
one can retrieve the artificial clusters from the catalogue 
using the exact same cluster selection procedure as used 
for the data. This procedure is beyond the scope of this 
paper, since we mainly want to demonstrate that there 
is a degeneracy between incompleteness effects and the 
mass independent disruption model of CFW06. 

An important consequence of our result is that it con- 
firms that cluster life-times are strongly dependent on 
the environ ment in which t hey evolve, which was already 
suggested bv lHodgel (|1987l ) based on the difference in the 
age distribution for clu sters in the SMC a nd the solar 
neighbourhood (see also lLamers et aUl2005l for a compar- 
ison of cluster life-times in f our galaxies). In contrast to 
this, IWhitmore et all (|2007t l use the result of CFW06 to 
show that the cluster dN/dt distributions in the Anten- 
nae galaxies and in the SMC are very similar. Although 
it would simplify things if all clusters evolve similarly 
the first Gyr, it is rather counter-intuitive since then the 
evolution would not depend at all on the strength of the 
tidal field and the number of giant molecular clouds. Our 
results support a scenario where the life-time of clusters 
that have survived the gas removal phase is determined 
mainly by environmental factors. 
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